Optimal ferrofluids for magnetic cooling devices

Superior passive cooling technologies are urgently required to tackle device overheating, consequent performance degradation, and service life reduction. Magnetic cooling, governed by the thermomagnetic convection of a ferrofluid, is a promising emerging passive heat transfer technology to meet these challenges. Hence, we studied the performance metrics, non-dimensional parameters, and thermomagnetic cooling performance of various ferrite and metal-based ferrofluids. The magnetic pressure, friction factor, power transfer, and exergy loss were determined to predict the performance of such cooling devices. We also investigated the significance of the magnetic properties of the nanoparticles used in the ferrofluid on cooling performance. γ-Fe2O3, Fe3O4, and CoFe2O4 nanoparticles exhibited superior cooling performance among ferrite-based ferrofluids. FeCo nanoparticles had the best cooling performance for the case of metallic ferrofluids. The saturation magnetization of the magnetic nanoparticles is found to be a significant parameter to enhance heat transfer and heat load cooling. These results can be used to select the optimum magnetic nanoparticle-based ferrofluid for a specific magnetic cooling device application.

www.nature.com/scientificreports/ fluidic and magnetic behavior, which provides the advantage of controlling ferrofluid motion by an external magnetic field 23,24 . A stable ferrofluid usually implies that the size of the suspended MNP should be less than 15 nm in diameter 25 . Under the influence of a magnetic field, the MNP orient along the magnetic field direction, and the ferrofluid behaves like a liquid magnet. The extent of ferrofluid magnetization depends on the type (metallic, alloy, or ferrite) and concentration of MNP constituents in the ferrofluid, the strength and orientation of the applied magnetic field, and the temperature of the ferrofluid. Ferrofluid motion can be controlled by changing the magnitude and direction of the external magnetic field. Since the ferrofluid can be remotely manipulated in a wireless fashion, it has a multitude of applications in targeted drug delivery 26,27 , hyperthermia [28][29][30] , droplet microfluidics 31,32 , particle synthesis 33 , energy harvesting 5,34 and waste heat recovery [35][36][37] . Ferrofluid flow, under the combined effect of both temperature and magnetic field gradients, can be controlled by an external magnetic field. This phenomenon is called thermomagnetic convection 23 and is used to develop passive magnetic heat transfer systems. Moreover, magnetic cooling devices based on the thermomagnetic effect are noise-free, vibration-free, require no or low maintenance, and are self-pumping and self-regulating. The use of ferrofluids as a coolant results in a thermal conductivity enhancement of ~ 300% under an external magnetic field 38,39 .
In a magnetic cooling device, the thermal gradient arising from the difference between the heat load and heat sink temperatures results in a magnetization gradient in the ferrofluid column near the heat load region. As a result, the ferrofluid at the heat load possesses low or zero magnetization, whereas the ferrofluid near the heat sink possesses higher magnetization. As a result, the differentially heated ferrofluid starts to flow around the closed loop in the presence of an external magnetic field due to the non-uniform magnetic volume force. A schematic of such a device is presented in Fig. 1.
However, the passive nature of magnetic cooling devices depends on the type of magnetic field used.
A. AC magnetic field such as an electromagnet requires external energy input. Hence, any magnetic cooling device using an AC magnetic field is a passive cooling system with low external energy. Magnetic cooling devices of various length scales viz., microscale (mm length-scale, mW power), small-scale (cm length-scale, W power), and large-scale (m length-scale, kW power), were developed and investigated for their thermomagnetic cooling performance. Earlier investigations mainly examined the cooling of microscale devices, e.g., for MEMS and LOC applications [40][41][42][43] . Small scale devices were researched to cool power electronics and small mechanical systems. Lian et al. 44 performed particle image velocimetry (PIV) measurements to determine ferrofluid velocity for several heat loads and flow channel geometries. Li et al. 45 developed a rectangular cross-section magnetic cooling device to examine the effect of external magnetic field strength and heat load power on the ferrofluid flow and heat transfer characteristics.
Recent studies focused on magnetic cooling devices with meter size flow channel length and larger channel diameter for long-distance and high-power waste heat transfer. Pattanaik et al. 46 performed kW level cooling of heat load by a multi-torus magnetic cooling device. The device was with a 1.8 m long effective channel, 25 mm channel diameter, and used a commercial oil-based Fe 3 O 4 ferrofluid as a working fluid. Khairul et al. 17 investigated magnetic field effects on thermomagnetic convection for laminar, turbulent flow conditions for a waterbased magnetite ferrofluid. Yamaguchi et al. 47 used a temperature-sensitive binary magnetic fluid to develop a 5 m perimeter magnetic cooling system. Goharkhah et al. 21 modeled a ferrofluid's forced thermomagnetic heat transfer under the influence of eight AC magnetic field line source dipoles. Fadaei et al. 48 numerically investigated the effect of magnetic field rotation on cooling performance. Hybrid ferrofluids were also investigated to improve the thermomagnetic cooling performance. Shahsavar et al. 49 experimentally investigated the heat transfer performance of CNT-loaded ferrofluid in the presence of constant as well as alternating magnetic fields. Sadeghinezhad et al. 50 experimentally studied the improvement in convective heat transfer using a graphene-loaded ferrofluid under an external magnetic field.
There are a few reports focused on the thermomagnetic cooling performance of ferrites-and metallic-based ferrofluids. Love et al. 40 designed and investigated a magnetocaloric pump for lab-on-chip applications. They synthesized oil-based magnetite and water and oil-based Mn-Zn ferrite ferrofluids and compared their cooling in a magnetic cooling system. Chaudhary et al. investigated the magnetic cooling performance of a water-based Mn-Zn ferrite ferrofluid using both experiments and simulations. They investigated the effect of heat load temperature, the concentration of magnetic nanoparticles in the ferrofluid, and saturation magnetization of the magnetic nanoparticles on the cooling performance. A similar research work using Mn-Zn ferrite was carried out by Phor et al. 51 . Sharma et al. 37 synthesized Fe-Cr-Al magnetocaloric ferrofluids and showed the effect of chromium content in Fe-Cr-Al, magnetic nanoparticle volume fraction, and heat load power on magnetic cooling.
These reports investigated the cooling performance as a function of device properties, heat load parameters, magnetic field strength, position, and orientation. Most of the literature described the use of ferrite-based commercial ferrofluids. For enhanced heat load cooling, Aursand et al. 52 numerically investigated cooling based on their optimization results for solenoid geometry, power consumption in the solenoid, thermomagnetic pumping, and MNP size.
Despite its apparent critical importance in cooling, there is surprisingly no formal ranking of the cooling performance of various ferrite-or metallic-MNP ferrofluids.
Hence, the scope of the present study is to investigate the cooling performance of the magnetic cooling device for a range of values of the magnetic and thermophysical properties of the ferrofluid. For the first time, we modeled the performance metrics, non-dimensional parameters, and heat load cooling performance to rank ferrite and metallic ferrofluids by their cooling performance. The magnetic pressure term was derived and found to be a function of saturation magnetization, magnetic nanoparticle volume fraction, temperature-dependent initial magnetic susceptibility, and magnetic field strength. The friction factor, power transferred from the heat load to the heat sink, and the exergy loss were derived as a function of magnetic, thermal, and thermophysical parameters. Numerical simulations were performed to examine the effect of magnetic parameters of the nanoparticles viz., bulk saturation magnetization, Curie temperature, pyromagnetic coefficient, and initial magnetic susceptibility on the cooling performance. Several ferrite and metallic ferrofluids were included in this study. www.nature.com/scientificreports/ It was found that ferrite-based ferrofluids of γ-Fe 2 O 3 , Fe 3 O 4 , and CoFe 2 O 4 exhibited superior cooling performance. For the case of metallic ferrofluids, FeCo based ferrofluid exhibited the best cooling performance, followed by Fe and FeNi ferrofluids. The bulk saturation magnetization of the MNP was found to be a significant parameter that enhanced the heat transfer and heat load cooling. Higher thermal conductivity and lower viscosity give rise to higher cooling performance. The Nusselt number, magnetic Rayleigh number, and Peclet number increased significantly for higher bulk saturation magnetization. Our ranking of various ferrofluids can be used to select the optimum ferrofluid for enhanced heat load cooling for a specific application.
Analytical and modeling methods "Numerical model" and "Governing equations for numerical model" sections describe the basics of the numerical model to investigate the thermomagnetic convection of various ferrofluids and the governing equations of the same. "Analytical methods for the derivation of the figure of merit" section presents a simple analytical derivation of the figure of merit using 1D governing equations for the thermomagnetic flow around a closed loop. "Nondimensional parameters and scaling" section briefs the non-dimensional parameters required to understand heat transfer in magnetic cooling devices. "Boundary condition of the numerical model" and "Numerical verification and validation" sections list all the boundary conditions and validation of the numerical model respectively, which are required to investigate the magnetic cooling by ferrite-and metallic-based ferrofluids. "Parameters considered" section includes all the parameters considered to simulate the thermomagnetic convection effect in this research work.
Numerical model. COMSOL Multiphysics-based 2D numerical simulations investigated the magnetic cooling performance of various ferrite-and metallic-based ferrofluids for a racetrack device. The ferrofluids were treated as a single-phase continuum as per Rosensweig's continuum approach 25 . Thermal, fluidic, and magnetic properties of the ferrofluids were provided as inputs to the numerical model. Three physics modules were employed, viz., magnetic field, heat transfer, and single-phase flow in the laminar regime. The magnetic volume force was added as an additional volume force term in the momentum conservation equation due to the thermomagnetic convection. Both magnetic volume force and the magnetic field vary non-linearly with distance (r) from the magnet 53 . To avoid any error that might arise from the discretization of magnetic volume force, we defined the corresponding components at the wall of the tube. Dense meshing was used along the boundary to avoid any errors due to the thin boundary layer. The magnetic field distribution of a suitable NdFeB magnet was simulated and compared to the experimentally obtained magnetic field strength values. The boundary conditions were employed as per the experimental design. The 2D numerical model was validated to obtain simulation accuracy.
Governing equations for numerical model. The following governing equations are used to model the thermomagnetic convection effect. The magnetic field distribution around the NdFeB permanent magnet was given by Maxwell's equations 54-56 : where H is the magnetic field, V m is the magnetic scalar potential, B is the magnetic flux density. μ 0 and χ m are the magnetic permeability of free space and magnetic susceptibility of the ferrofluid, respectively.
The non-uniform magnetic susceptibility of the ferrofluid as a function of temperature and the magnetic field is given by the Langevin eqaution 33,55,56 : where C 0 is the volume fraction of the magnetic nanoparticles in the ferrofluid, M Bulk s is the bulk saturation magnetization of the ferrofluid, and γ is the Langevin parameter.
L is the Langevin function, which takes the form: The Langevin parameter γ is a function of temperature and is defined as: where χ i is the initial susceptibility of the magnetic nanoparticles over a range of temperature values below the Curie temperature.
(1) www.nature.com/scientificreports/ The transient form of the continuity equation, the Navier-Stokes equation, and the heat transfer equation describe mass, momentum, and energy conservation 46 . These transient forms of the mass, momentum, and energy conservation equations are used in the model to simulate the thermomagnetic convection effect.
where T, v, and p denote temperature, velocity, and pressure of the ferrofluid, respectively. F m is the non-uniform magnetic volume force acting on the ferrofluid. The ferrofluid's thermophysical properties are dynamic viscosity ( η ), density (ρ), specific heat (C p ), and thermal conductivity (κ).
The magnetic volume force (F m ) experienced by the ferrofluid in the presence of an applied magnetic field is given by Pattanaik et al. 46 : where, χ mnp represent the volume magnetic susceptibility of the magnetic nanoparticles. χ m is the magnetic susceptibility of the ferrofluid.
Analytical methods for the derivation of the figure of merit. The analytical derivation presented in this section is based on a simple one-dimensional (1-D) approach (Fig. 2). The axis along the middle of the flow loop (l, flow loop length) is taken as the 1-D coordinate to derive the ferrofluid temperature, magnetic pressure (Δp m ), and the friction factor (f). The schematic of the 1-D flow channel loop that is used to derive the temperature of the ferrofluids at various regions of the flow loop, the magnetic pressure, and the friction factor is shown in Fig. 2a, b. The middle of the heat load section is taken as the origin of l. The l value at various points along the flow loop is labeled in the figure. The magnetic cooling device has 4 different regions: (a) heat load region, (b) outflow channel, (c) heat sink region, and (d) inflow channel, as shown in Fig. 2b. While developing the reduced 1-D governing equations, the cross-sectional average is considered for the thermomagnetic convection. The following assumptions were taken into consideration while deriving for the temperature, magnetic pressure, and friction factor.
• The cross-sectional average temperature of the ferrofluid is similar to the bulk temperature of the ferrofluid along the axial flow loop length (l). • Axial conduction of heat is negligible compared to convection-based heat transfer. Here, α is the thermal diffusivity of the ferrofluid.
The energy conservation equations for the heat load section, heat sink section, and the inflow and outflow channels can be written as, www.nature.com/scientificreports/ Here, P is the applied heat load power, and h HS is the heat transfer coefficient at the heat sink. The length of the heat load and the heat sink is denoted as l HL and l HS , respectively. T f and T w denote the temperature of the bulk ferrofluid and the flow channel wall, respectively. The right-hand side of the equation is zero because of the www.nature.com/scientificreports/ negligible temperature difference between the ferrofluid and the flow channel wall and the adiabatic nature of the flow channels. Integrating the above energy equations lead to the following temperature expressions, The above expressions provide the ferrofluid temperature at any point along the flow channel. Equations (19)- (22) are the expressions for evaluating the ferrofluid temperature as a function of the flow loop length (l) at the heat load region, along the outflow channel, heat sink region, and along the inflow channel, respectively. T in , T out , and T HS are temperatures of ferrofluid at the inlet, the outlet temperature, and the heat sink, respectively. l ′ represents the length of the inflow and outflow channels of the racetrack magnetic cooling device.
The 1D momentum equation (Eq. (14)) can be integrated around the racetrack flow channel to obtain the friction factor as a function of the resultant magnetic pressure due to the thermomagnetic convection effect. Integration of the steady-state momentum equation takes the form, The integration of the convective acceleration term (1st term in Eq. (23)) is zero owing to the continuity equation (Eq. (13)), under the assumption of the constant density of the ferrofluid across the flow channel. The pressure gradient term around the flow channel vanishes as the pressure difference along a closed loop is zero. Integration of the viscous drag term (3 rd term in Eq. (23)) leads to Eq. (24) 57 . where f is the friction loss associated with the magnetic cooling device. l d and d t are the total perimeter and internal diameter of the flow channel, respectively. Integrating the added magnetic volume force term in the Navier-Stokes equation gives rise to the resultant magnetic pressure (Δp m ) due to the thermomagnetic convection effect.
H max and H min correspond to the maximum and the minimum field experienced by the ferrofluid in the vicinity of the heat load region. l max and l min are the corresponding positions along the flow channel with respect to the maximum and minimum applied magnetic field, respectively. In the limit of zero magnetic fields, Hence, the magnetic pressure term is negligible at low applied magnetic fields and vanishes at zero magnetic field. The magnetic pressure term is directly proportional to the ferrofluid's bulk saturation magnetization and magnetic nanoparticle volume fraction. The temperature dependence arises from the Langevin parameter (γ).
Following the substitution of Eqs. (24) and (26) in Eq. (23), we obtain the friction factor (f) 17 , Exergy refers to the maximum useful work a system can perform when it is brought into equilibrium with its surrounding 58,59 . Exergy loss refers to the irreversible lowering of the maximum available work outside the control volume 17 . For thermomagnetic convection based magnetic cooling device, under the assumption of www.nature.com/scientificreports/ zero energy exchange between the system and the surrounding, the amount of exergy loss at the steady-state of a ferrofluid is given by Khairul et al. 17 , Here, T amb refers to the ambient temperature. • Use of soft magnetic nanomaterials with high magnetic permeability and extremely low coercivity for faster response to the external magnetic field • High pyromagnetic coefficient near the heat load temperature for superior heat load cooling • High saturation magnetization • Low magneto-crystalline anisotropy of magnetic nanoparticles 53 • Low viscosity of the ferrofluid The first 4 refers to the properties of the magnetic nanoparticles. We have chosen a range of ferrite-and metallic-based nanoparticles based on the criteria. Table 1 summarizes the parameters, their units, notation, and experimental sets considered in this work.

Non
(29) P ′′ = T amb ṁC p ln T out T in Table 1. Parameters, their notation, range, and the measurement sets.

Results and discussions
Experimental and simulated heat load cooling for magnetite ferrofluid. Before performing a numerical investigation for various ferrofluids, the simulation results were validated against the experimental cooling results obtained from a copper racetrack-shaped magnetic cooling device (Fig. 1) with a perimeter of 130 cm and an internal tube diameter of 0.81 cm. T-type thermocouples were mounted on the surface of the heat load at various locations, as shown in Fig. 1. Data acquisition was performed by a lab-view setup and a data logger (Pico technology, model: TC-08, channels: 8).
The experimentally obtained heat load temperature profile was compared with the simulated profile for a heat flux value of 3.49 kW/m 2 and an initial heat load temperature of 197 °C for a copper-based racetrack magnetic cooling device containing commercial Fe 3 O 4 ferrofluid. The experimental and simulated temperature profiles are in good agreement (Fig. 3). This simulation was employed to study the effect of ferrofluidic thermophysical and magnetic properties on cooling performance.

Parametric investigation of heat load cooling.
A preliminary study to assess the significance of various magnetic parameters of the ferrofluid viz., initial magnetic susceptibility (χ i ), pyromagnetic coefficient (dχ/ dT), bulk saturation magnetization ( M bulk s ), and the Curie temperature (Tc) on the heat load cooling was conducted. A reference susceptibility versus temperature graph was constructed. The experimental values of the magnetite (Fe 3 O 4 ) nanoparticles χ′ (T) were used as a guide 60 . A Tc of ~ 200 °C was selected to be close to the range of heat load temperatures investigated (Fig. 5a). The results are described in the following subsections.
Effect of bulk saturation magnetization on heat load cooling. For the reference susceptibility curve (Fig. 4a) and an initial heat load temperature of 100 °C, the effect of the M bulk s magnetic nanoparticles on heat load cooling performance and average ferrofluid velocity were investigated, as shown in Fig. 4b, c. ΔT increased significantly from 2 °C to 52 °C with an increase of Ms from 50 to 500 kA/m. It followed a polynomial behavior. The average ferrofluid velocity also increased significantly in a linear manner as a function of bulk Ms. The higher value of bulk saturation magnetization results in larger magnetic pressure (Eq. (26)), which gives rise to higher ferrofluid velocity, and hence better cooling.
Effect of change in initial susceptibility on heat load cooling. The effect of initial magnetic susceptibility (χ i ) on the heat load cooling was investigated by assuming two different χ i profiles with the same pyromagnetic coefficient and same Curie temperature of 200 °C (Fig. 5a,b). These susceptibility profiles were provided as inputs to the numerical model. Subsequently, the effects of heat flux and the bulk saturation magnetization on the cooling performance were investigated for a heat load temperature of 200 °C. Figure 5c, d shows the effect of heat flux on the heat load cooling and the average ferrofluid velocity for these two susceptibility profiles for a bulk saturation magnetization value of 446 kA/m. The heat load temperature was 200 °C. Both the heat load cooling and average ferrofluid velocity increased with larger heat flux values. However, for temperatures higher than the Tc of the ferrofluid, cooling and the velocity reduced sharply. This ΔT reduction is due to the drastic reduction of the thermomagnetic force in the paramagnetic regime of the ferrofluid. The increased cooling for higher heat flux values (below Tc) reveals the self-regulating nature of such magnetic cooling devices.
We also investigated the effect of the ratio of initial magnetic susceptibilities to the reference susceptibility profile on the cooling performance (Fig. 6). Various χ i vs. temperature profiles were considered, as shown in Fig. 6a. Both the heat load cooling (Fig. 6b) and the ferrofluid velocity (Fig. 6c) showed a slight increase of 10.6% www.nature.com/scientificreports/ and 29%, respectively, with a 700% increase in n value (χ″/ χ′, χ′ being the reference χ curve). Hence, initial magnetic susceptibility did not significantly affect cooling.
Effect of change in pyromagnetic coefficient on heat load cooling. The effect of the pyromagnetic coefficient on the cooling performance was investigated by several virtual susceptibilities vs. temperature curves with increasing pyromagnetic coefficient ratio (Fig. 7a). The corresponding pyromagnetic coefficient vs. temperature curves are shown in Fig. 7b. The dependence of the heat load cooling magnitude and the average ferrofluid velocity on the pyromagnetic coefficient ratio is plotted in Fig. 7c, d. Both the curves exhibited an increasing trend and showed similar behavior. However, the cooling performance tends to saturate for a higher ratio of the pyromagnetic coefficient. The heat load cooling magnitude increased by only 4 °C for a 200% increase in the pyromagnetic coefficient. Hence, the pyromagnetic coefficient is also not a significant parameter to enhance the cooling.
Effect of Curie temperature on heat load cooling. To study the effect of Curie temperature (Tc) on magnetic cooling of the heat load, we considered several virtual χ i vs. temperature curves with the same values of M s and initial χ i (Fig. 8a). The effects of Tc on heat load cooling and the ferrofluid velocity profile were studied for a heat load temperature of 60 °C and for a range of bulk saturation magnetization values (Fig. 8b,c). Both of these parameters exhibited an increase in the cooling for values of Tc closer to the heat load temperature. In the present case, a Tc value of 75 °C is closest to the heat load temperature; hence, it exhibited maximum cooling performance (Fig. 8b) as well as the highest velocity (Fig. 8c). These results can be attributed to the high pyromagnetic coefficient (dχ/dt| T=60 °C ) when the Tc is closer to the heat load temperature. This result is consistent with the effect of the pyromagnetic coefficient (Fig. 7c,d).
From these results, it is clear that M bulk s is one of the most significant parameters to improve cooling. For constant M bulk s , the dχ/dt and Tc affect the cooling. On the other hand, the χ i has a negligible effect on the heat load cooling performance. The effects of these parameters can also be understood from the derived expression of magnetic pressure (Eq. (26)). were performed using the initial susceptibility vs. temperature curve to rank the ferrite and metallic nanoparticle-based ferrofluids with respect to their thermomagnetic cooling performance. The bulk saturation magnetization of the nanoparticles was provided as input to the model. Figure 9a, b shows the heat load cooling and the average ferrofluid velocity of various ferrite-based ferrofluids as a function of the bulk saturation magnetization of corresponding magnetic nanoparticles. Figure 9c, d is the heat load cooling and the average ferrofluid velocity curves for metallic nanoparticle-based ferrofluids. It shows an increase in heat load cooling and ferrofluid velocity with increasing bulk saturation magnetization for both cases. The cooling and ferrofluid velocity significantly enhanced ferrofluids having higher bulk Ms at high heat flux and heat load temperature. Certain anomalies were observed for the ΔT and the velocity curves (yellow shaded area) in the case of ferrite-based ferrofluids. These anomalies can be attributed both to the pyromagnetic coefficient and the Curie temperature values. At higher temperatures, the pyromagnetic coefficient is larger. Hence, the anomalies are significant at larger heat load temperatures at higher heat flux values. A decrease (1 st and 3 rd yellow shaded regions) or increase (2 nd and 3 rd yellow shaded regions) in cooling and the velocity corresponds to a change in the dχ/ dT value. A slight anomaly is observed in the lower bulk Ms region for the metallic ferrofluids case. However, due to the high bulk saturation magnetization of metallic ferrofluids, ΔT and average velocity are significantly larger in comparison to the ferrite ferrofluids. Tables 2 and 3 summarize the cooling performance for various ferrite and metallic ferrofluids, respectively.
Hence, at lower heat load temperatures, ferrite ferrofluids having moderate magnetization values exhibit similar cooling performance as their metallic counterparts. On the other hand, at higher heat load temperatures, metallic ferrofluids are more attractive. However, it is usually challenging to synthesize stable metallic magnetic nanoparticle-based ferrofluids due to the oxidation of the metal and the high density of the metal nanoparticles. These issues may lead to particle agglomeration, settling, and a reduction in magnetization of metallic ferrofluids.
Non-dimensional parameters and figure of merit of ferrofluids. Transferred power and exergy loss versus bulk saturation magnetization. Figure 10a shows a plot of waste heat power vs. bulk saturation mag- www.nature.com/scientificreports/ netization of the magnetic nanoparticles. Waste heat power is the power transported from the heat load to the heat sink by a magnetic cooling device due to the thermomagnetic convection of the ferrofluid. More waste heat power can be removed using highly magnetic ferrofluids. This is due to the higher volume flow rate of ferrofluids having higher magnetization for given magnetic field strength. Figure 10b shows the exergy loss curve versus the bulk saturation magnetization. The exergy loss also tends to increase with bulk Ms. Exergy loss represents higher randomness. Hence, the entropy maximizes at higher exergy loss. Both the power and the exergy loss are directly proportional to the mass and volume flow rate of the ferrofluid, which is larger for higher bulk Ms due to the greater thermomagnetic volume force experienced by the ferrofluid. The ferrofluid forms vortices near the heat load region due to the combined effect of the thermal gradient and the magnetization gradient 46 . These vortices become prominent for a high magnetization ferrofluid. Hence, the increase in entropy can also be attributed to the formation of strong vortices when the ferrofluid possesses higher magnetization.

Non-dimensional parameters versus bulk saturation magnetization.
We examined the behavior of several nondimensional parameters as a function of the bulk saturation magnetization of magnetic nanoparticles to understand the thermomagnetic convection-based heat transfer of a magnetic cooling device.
The nature of heat transfer at the heat load region was investigated by calculating the average Nusselt number (Nu avg ) using Eq. (5), as given in the Supplementary file. Nu avg increases for larger bulk M s (Fig. 11a), implying greater thermomagnetic convection heat transfer from the heat load. This observation is also consistent with the increase in the magnitude of cooling (ΔT) of the heat load for higher bulk M s (Fig. 9). Nu avg is below 10 for bulk M s values below 250 kA/m, suggesting that the ferrofluid flow is laminar. However, for bulk M s values above 250 kA/m, the ferrofluid flow can be considered as a transition between the laminar and turbulent flow regime. The red curve in Fig. 11a suggests a linear increase in the Nu avg with bulk M s . Hence, the passive thermomagnetic convection heat transfer increases for ferrofluids having higher saturation magnetization.
The Ra m increases for larger values of bulk Ms (Fig. 11b), which is consistent with Eq. 6 of the Supplementary file. It signifies a stronger dominance of the thermomagnetic force vs. the viscous force for magnetic cooling systems. The Ra m of 10 9 and above signifies the stronger turbulence effect due to the formation of ferrofluid vortices. www.nature.com/scientificreports/ Figure 11c shows the Peclet number (Pe) as a function of bulk M s . The linear increase in Pe with increasing bulk M s suggests greater advection heat transfer compared to diffusion-based heat transfer. Thus, the convective motion of the bulk fluid carries a significant amount of waste heat away from the heat load region compared to conduction heat transfer through the nanoparticle chains. Figure 11d shows the effect of Stanton number (St) as a function of bulk M s . The non-linear reduction in St with an increase in bulk M s signifies a reduction in the thermal capacity of the ferrofluid. Conversely, the thermal capacity tends to increase for a higher volume flow rate of the ferrofluid. Figure 12 shows the plot of the average Nusselt number against the magnetic Rayleigh number. Nusselt number increases with an increase in Ra m value. As expected, the thermomagnetic force is higher for higher Ra m values, giving rise to a stronger thermomagnetic convection effect, resulting in higher average Nu. However, the average Nu shows two distinct trends. At lower values of Ra m , Nu scales as Ra 0.89 m , whereas at higher Ra m values, Nu scales as Ra 0.62 m . This can be attributed to the viscosity of the ferrofluid at the heat load region. When the ferrofluids with higher saturation magnetization (e.g., metallic ferrofluids) are considered, the heat load is cooled significantly compared to lower magnetization ferrofluids. This heat load cooling results in lower ferrofluid temperature at the heat load region for high magnetization ferrofluids compared to low magnetization ferrofluids. Hence, the high magnetization ferrofluids exhibit larger effective viscosity, resulting in lower thermomagnetic convection. This interplay between the magnetization and the viscosity results in two different regimes.

Conclusions
The cooling performance of the thermomagnetic convection-based magnetic cooling device was numerically assessed with respect to the magnetic properties of various ferrite-and metallic/alloy-nanoparticle-based ferrofluids. For the first time, we ranked the ferrofluids based on their thermomagnetic cooling performance. The effect of bulk saturation magnetization, initial magnetic susceptibility, pyromagnetic coefficient, and the Curie temperature on the heat load cooling was evaluated. Several ferrites-and metallic/alloy-based ferrofluids were www.nature.com/scientificreports/ simulated to investigate their thermomagnetic cooling performance. Various performance metrics and nondimensional parameters were evaluated. The conclusions are as follows.
• The bulk saturation magnetization plays the most significant role in cooling the heat load. The initial magnetic susceptibility and pyromagnetic coefficient are less significant. • Maximum cooling performance is obtained when the Curie temperature is near the heat load temperature. • γ-Fe 2 O 3 , Fe 3 O 4, and CoFe 2 O 4 ferrofluids displayed better cooling performance than other ferrite-based ferrofluids. • For metallic/alloy-based ferrofluids, FeCo ferrofluid exhibited the best cooling performance, followed by Fe and FeNi ferrofluids. • Power transport capability and exergy loss increased for higher bulk saturation magnetization of magnetic nanoparticles. • The non-dimensional magnetic Rayleigh number, average Nusselt number, and Peclet number confirmed the enhancement in thermomagnetic convection-based passive heat transfer for strongly magnetic ferrofluids.     www.nature.com/scientificreports/